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Abstract: 

It is interesting to study the characteristics of the whole fam- 
ily of B c which contains two different heavy flavors. LHC and 
the proposed Z° factory provide an opportunity because a large 
database on the B c family will be achieved. B c and its excited 
states can be identified via their decay modes. As suggested by 
experimentalists, B*(ns) — > B c + 7 is not easy to be clearly mea- 
sured, instead, the trajectories of n + and n~ occurring in the decay 
of B c {ns) B c (ms) + 7r + 7r~ ( n > m ) can be unambiguously 
identified, thus the measurement seems easier and more reliable, 
therefore this mode is more favorable at early running stage of 
LHCb and the proposed Z° factory. In this work, we calculate the 
rate of B c (ns) — > B c (ms) + tt + tt~ in terms of the QCD multipole- 
expansion and the numerical results indicate that the experimental 
measurements with the luminosity of LHC and Z° factory are fea- 
sible. 

PACS numberfs): 13.25.-k 
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I. INTRODUCTION 



The heavy quarkonia such as charmonia and bottomia have been experimentally and 
theoretically explored for several decades already because large database on them is available. 
In comparison, the physics on B c has not been thoroughly studied yet. The reason is obvious 
that B c contains two different heavy flavors, so unlike quarkona, it is produced via more 
suppressed processes at e + e~ colliders. The earlier work [l| indicates that at the luminosity 
of regular e + e~ colliders (the LEP I and II), one cannot expect to observe B c production, 
i.e. its production rate is too small to be measured. Therefore, people should turn to hadron 

n n 

colliders. As predicted B c was observed at TEVATRON[2| a while ago. At the energy and 
luminosity of LHC, one may anticipate B c events thousand times more than at TEVATRON. 
Now a project of constructing a Z° factory is proposed which will provide sufficiently high 
luminosity at the Z° pole. Even though the production process e + e~ —> B c + B c where 
a pair of heavy quarks (cc or bb) emerges from a hard gluon emission, is suppressed, the 
high luminosity and the pole effect would greatly enhance the production rate, i.e. the 
enhancement compensates the suppression and enables the production of B c measurable. 

Since B c is made of two different heavy flavors, its decay characteristics would somehow 
distinct from the heavy quarkonia which contains a pair of heavy quarks of the same flavor. 
Namely, the two constituents in quarkonia annihilate into gluons which afterwards hadronize. 
Although this mode is OZI suppressed, it is a strong interaction-induced process and has 



B c family can only decay via weak 
3|]. It is interesting to investigate such 



a larger decay width. Nevertheless the ground state of 
interaction and its lifetime has been carefully studied 
mesons and we are not only interested in the ground state, but also its excited states. The 
lowest excited states would be the vector B* (Is) and pseudoscalar B c (2s) (0~) and the 
latter one is the first radially excited state of the family. A simple analysis estimates the 
splitting between B c (ls) and B*(ls) is about 50 to 80 MeV, so that B* c ->• B c + 7r° is 
forbidden by the energy-conservation and the hadronic decays of B* can only occur via weak 
processes which are obviously suppressed. Only possible transition is the radiative decay 
B* — > B c + 7, but at LHC, detection of a single photon from a messy background is very 
difficult. Instead, in the decay B c (ns) — > B c (ms) + 7r + 7r~ the daughter charged pions are 
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very easy to be identified. Therefore, our experimental colleagues strongly suggest us to 
investigate the channel B c (ns) — > B c (ms) + 7r + 7r~. Definitely, we can gain more information 
about the B c family. 

Since B c is composed of two heavy quarks, the relativistic effects may not be too serious, 
thus the potential model can be a good choice for determining the spectra of B c and its 
excited states. In analog to dealing with heavy charmonia and bottomia, we employ the 
Cornell potential to calculate the masses of B c (ns) and B*(ns) where the mass of the ground 
state B c (ls) is taken as input for fixing concerned model parameters. 

Following literature 4]-|8| we evaluate the decay rate of B c (ns) — > B c {ms) + tttt in terms 
of the QCD multipole-expansion. The picture is depicted as that the initial B c transits 
into a hybrid state beg where be stays in a color-octet, by emitting a gluon, and then the 
hybrid turns into B c (ls) by emitting the second gluon. The two gluons eventually hadronize 
into two pions. In the transitions B c (ns) — > B c (ms) + ti + ir, the momentum transfer 
is not large and the perturbative method does not apply. 
(QCDME) method suggested by Gottfried, Yan and Kuang 
meson emission process. In the picture of the QCD multipole expansion, the two emitted 
gluons are not energetic particles, but described by a chromo filed of TM or TE modesQ]. 
It is worth emphasizing again that the two gluons are not free gluons in the sense of the 
perturbative quantum field theory, but a field in the QCD multipole expansion. It is easy to 
understand that such transition is dominated by the El-El mode, while the Ml-Ml mode 
is suppressed for the heavy quarkonia case. 

The two chromo-El transitions are dealt with in the regular framework of multipole- 
expansion. The key point is how to properly obtain the mass-spectrum and wavefunction of 
the intermediate hybrid state. Isgur and Paton[l(| suggested to use a modified potential to 

3ut there are a few free parameters to be determined. Before, the 



?he QCD multipole expansion 
4-^] properly treats the light- 



describe the hybrid states 
parameter in the potential llj which Kuang and Yan employed, was fixed by assuming the 
^(4040) to be the hybrid 5j. Thanks to the achievements of BELLE, CLEO and BES, an 
abundant database on such two-pion emission processes has been available. In our previous 
work Ijj], we carefully discussed the cases of T(ns) — > T(ms) + 7r + 7r and tp(ns) — > t/j(ms) + 
Ti + Ti in terms of the QCD multipole-expansion, then by the method of minimizing \ 2 which 



is widely adopted in analyzing experimental data, we eventually fix the parameters in the 

Q fl and 



potential. For the hybrid state we used three differe nt p otential models 
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noticed that the potential form proposed by Allen eta/.[14j better coincides with the lattice 
result. According to the the potential parameters gained by fitting the spectra of heavy 
quarkonia T and ip families, one can estimate the parameters for the members of the B c 
family, by slightly varying the values of corresponding parameters (see below for details). 

In this work, we apply the QCD multipole expansion method 5[ and the potential form 
suggested by Allen et al. jl^] to calculate the transition rate of B c (ns) — > B c (ms) + ir-K where 
the potential parameters for the hybrids \bcg > (\bcg >) are the same as that we obtained 
for \bbg > and \ccg >. Since the spectra and wavefunctions of higher exited states of B c are 
even harder to be accurately derived, at this stage, we only concern the transitions from the 
radially excited states B c (3s) and B c (2s) into the lower states via emitting two pions. 

The derivation and numerical computations are straightforward and very similar to the 
procedures carried out in literature, therefore, below unless necessary for clarity, we omit 
some technical details. This work is organized as follows. After the introduction, we present 
the theoretical formulae for the transition B c (ns) — > B c (ms) + 7T7T, and then in sec. Ill, we list 
our numerical results along with all input parameters in tables, the last section is devoted 
to our discussion and conclusion. 

7T+ + tt-{2tt°) 




FIG. 1: The two-gluon emission diagram and the two gluons eventually hadronize into two pions. 
The intermediate state is a hybrid state \bcg >. 
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II. FORMULATION 



The transition width 



The theoretical framework about the QCD Multiploe Expansion method is well framed in 
Refs 5|-|8|, and all the concerned formulas are presented in those papers. Here we only copy 
a few formulas which are necessary for evaluating the widths of B c {ns) — > B c (ms) + 7r + 7r 
(n > m) in this work. For all details, the readers are suggested to refer to the original works 
and references therein. The general formula for the rate caused by double El transitions was 
given in Ref. 5| as 

' I 1 k 


s h J f 



- 8 lilf 5 JiJ ^G\C l \ 2 -\h\C^) |£(2/ + 1) ' 

'000 



+ (2k + 1)(2/ / + 1)(2J / + 1) J2( 2k + !)[! + (-1; 



V 



f! 



I |2 



H\C,\ 2 



£(2/ + 1) 



f h 1 l \ 


( 1 1 h\ 


h I 1 


. ) 


K o o oj ' 


k 1 k If [ 



)fi 



I |2 
i,f\ ' 



where and Jj(/) are orbital and total angular momenta of the initial( final) heavy quarko- 
nia and the total spin of the initial and final states respectively. / is the angular momentum 
of the color-octet intermediate state. |Ci| 2 and | C 2 1 2 are constants to be determined which 
come from the hadronization of gluons into pions. G and H are the phase-space integrals 
whose concrete forms were given in 5|, |6|. Obviously the first term corresponds to an S-wave 
and the second term to a mixing of S and D-waves. 

In terms of Eq.([T]) the transition rate of a pseudoscalar meson into another pseudoscalar 
meson with a two-pion emission can be written as 



T(n l 1 S Q ^n f l S ) = \C l \ 2 G\fl f 



l |2 



(2) 



which is similar to that between two vector-quarkonia[5|, |6(. ni,rif are the principal quantum 
numbers of initial and final states and fh is the overlapping integration over the concerned 
hadronic wave functions, 

J R f (r)r p fR* Kl (r)r 2 dr J Rl^r'y^Riir'y^dr' 

K 



fi,f - Yj 
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5 
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Kl 



(3) 



where Pi, Pf are the indices related to the multipole radiation, for the El radiation Pj, -P/=l 
and 1 = 1. Ri,Rf and Rki are the radial wave functions of the initial and final states, Mi is 
the mass of initial quarkonium and Eki is the energy eigenvalue of the intermediate hybrid 
state. The sum over the principle number of the intermediate hybrid state K is truncated at 
K = 7 because the contributions from higher excited states are too small and can be safely 
neglected. 

Because strong interaction is blind to flavor and electrical charge, the whole scenario 
should be applicable for the B c cases. 



B. The phenomenological potential 



In this study we employ generalized Cornell potential 15] which includes a spin-related 



term 



16] for the initial and final states as 



V(r) = + br + V s {r) + V , 

r 



(4) 



where k = 4ci3 ^ and the coupling a s {r) can be treated as a phenomenological constant while 
calculating the spectra of quarkonia. Thus, for the phenomenological application, one does 
not need to consider their QCD running. The spin-related term V s is, 

8ttk 



r, 



3m 2 Q 



(5) 



with 
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(6) 



and Vq is the zero-point energy 



12| 



The intermediate state between the two gluon-emissions is a hybrid state, namely the QQ' 
resides in a color-octet. It was indicated that one still can use an effective potential to describe 



the color-octet QQ' state |1Q| . 
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HQ 



In literature, there are four different effective 



potential forms which are respectively suggested by the authors of Refs. |10l 



11 



13 



14j. 



In our earlier work 12[ we emp 
\QQ'g > (Q stands as b or c in 



oyed the three models of them to study the hybrid state 



121 ] ) whereas, Yan and Kuang 



5( used the potential given 
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by Buchmiiller and Tye 



We find that the potential form suggested by Allen etal. 
coincides with the lattice result better than the others, thus we will use that potential in 
this work. 



The corresponding potential form is 



Vi(r) = - + v /(6V) 2 + 27r6' + y '. 



(7) 



Because the authors of Refs. 
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did not consider the spin-related term, we have 



modified the potential by adding the spin-related term V s which has the same form given 
above in Eq.(j5]), then the potential becomes: 

V(r) = V { + V.. (8) 



With this modification, we can investigate the spin-splitting effects. 



III. NUMERICAL RESULTS 

At first we need to fix the parameters in the potentials. In our earlier workfl^. we 
re-fitted the spectra of the quarkonia to obtain the corresponding potential parameters in 
Eq. fl4]) . The values of the parameters are listed in Table [H 



TABLE I: The potential parameters for cc and bb 





K 


6(GeV 2 ) 


m(GeV) 


cr(GeV 2 ) 


^o(GeV) 


cc 


0.67 


0.16 


1.78 


1.6 


-0.6 


bb 


0.53 


0.16 


5.13 


1.7 


-0.6 



One can see that the parameters b and Vq in Eq. (j4j) are the same for bb and cc, so that we 
suppose they are unchanged for be. Since the difference between the values of the parameter 
u for bb and cc is small, it is plausible to choose a = 1.65 for be (be). In the calculation for 
quarkonia nig is the mass of the constituent quark (b or c), instead, for be one should use 
the reduced mass. Here we set m c = 1.78 GeV and nib = 5.13 GeV for concrete numerical 
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TABLE II: our prediction on mass of some be mesons 





Is 


2s 


3s 


O-(GeV) 


6.276 


6.880 


7.254 


1- (GeV) 


6.356 


6.908 


7.274 



computations. Fitting the mass of B c (0~) 6.276 GeV (l7j, we obtain k = 0.58. Then with 
the given potential we predict the masses of a few other be states listed in the following table. 

For the hybrid potential (J7j), the strategy in our earlier work [3| is that we use the 
minimal y 2 method which is widely adopted in analysis of experimental data, to determine 
the potential parameters which are listed in Table HTT1 One can see that b' and Vq in Eq.(J7]) are 
the same but k' is different for bbg and ccg. Because no direct measurements on transitions 
(ms) + 7r + 7r have ever been conducted so far, we cannot determine the 
parameter k' for beg in terms of available data as we did for the heavy quarkonia. However 
one can expect that k' should fall in the region between the values for the bbg and ccg systems, 
thus we will vary this value slightly within the range to study a possible dependence of the 
numerical results on k'. The dependence is shown in Table HVl (see below). 

TABLE III: potential parameters for ccg and bbg 





k' 


b' (GeV 2 ) 


m(GeV) 


Vo'(GeV) 


ccg 


0.54 


0.24 


1.78 


-0.8 


bbg 


0.40 


0.24 


5.13 


-0.8 



Then there is still a free parameter |C*i| 2 in the decay rate Eq.(j2J). It is noted that C\ 
is a factor related to the hadronization of gluons into two pions, so should be universal for 
ip, T and B c meson decays. Because C\ is fully determined by the non-perturbative QCD 
effects, it cannot be derived from an underlying principle so far. Thus in Ref.fl^ we fixed 
| Ci | 2 = 182.12 x 10" 6 in terms of the well measured decay width T(ip(2S) — > J/ip + tx + tt). 

With these potential parameters, we solve the Schrodinger equations to obtain wave func- 
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tions and masses of the initial, final and intermediate states which appear in the overlapping 
integration fh. Then we can compute the corresponding widths of the concerned modes. As 
indicated above, in our calculation we vary k' in a small region from 0.40 (n'(ccg)) to 0.54 
(K'(bbg)). The corresponding results can be seen in Tabl dlVl For these numerical results, we 
observe that T(B C (3S) — > B c (2S))irir and T(B C (2S) — > B c (lS))inr are not very sensitive to 
the value of k', but T(B C (3S) — > 5 c (lS'))7r7r is. 



TABLE IV: prediction at a = 1.65 (widths in units of KeV) 





K 








0.40 


0.42 


0.44 


0.46 


0.48 


0.50 


0.52 


0.54 


M 6c - g (GeV) 


[l 


= 


K = 


1) 


7.629 


7.630 


7.631 


7.632 


7.633 


7.634 


7.635 


7.636 


M b5g (GeV) 


[I 


= 1 


K = 


1) 


7.800 


7.800 


7.801 


7.802 


7.803 


7.804 


7.805 


7.806 


M bSg (GeV) 


[I 


= 1 


K = 


2) 


8.118 


8.119 


8.120 


8.121 


8.121 


8.122 


8.123 


8.123 


M b5g (GeV) 


[l 


= 1 


K = 


3) 


8.417 


8.418 


8.418 


8.419 


8.420 


8.420 


8.421 


8.422 


M b5g (GeV) 


[I 


= 1 


K = 


4) 


8.700 


8.701 


8.701 


8.702 


8.703 


8.703 


8.704 


8.704 


M b5g (GeV) 


[I 


= 1 


K = 


5) 


8.970 


8.971 


8.971 


8.972 


8.972 


8.9735 


8.973 


8.974 


M b5g (GeV) 


[I 


= 1 


K = 


6) 


9.230 


9.230 


9.230 


9.231 


9.231 


9.232 


9.232 


9.233 


M b5g (GeV) 


[I 


= 1 


K = 


7) 


9.479 


9.480 


9.480 


9.481 


9.481 


9.481 


9.482 


9.483 


T(B C (3S) - 




B c (2S))7nr 


11.01 


11.07 


11.14 


10.90 


10.95 


10.99 


11.06 


11.11 


T(B C (3S) - 




B C (IS))TT7T 


4.91 


3.95 


2.98 


5.64 


4.85 


4.01 


3.05 


2.38 


T(B C (2S) - 




B C (IS))TT7T 


64.13 


64.01 


63.90 


63.71 


63.60 


63.48 


63.37 


63.25 



IV. OUR CONCLUSION AND DISCUSSION 

Study on B c , B* mesons and their radial and angular excited states is important because 
they are the last heavy mesons and are composed of two different heavy flavors. Because of 
their special structures, a thorough study on the production and decay processes where B c 
and its excited states are involved may shed more light on the fundamental interactions, espe- 
cially the non-perturbative QCD, moreover, may provide some hints to new physics beyond 
the standard model. Therefore, this field attracts attentions of theorists and experimental- 
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ists of high energy physics. The main obstacle for the study is that the production rate of B c 
is small as Chang and his collaborators indicated However, as LHC begins running, the 
high luminosity would provide sufficiently large data sample, moreover, a proposed Z-factory 
with a luminosity much higher than the LEP-I, would offer a clean environment for the B c 
research. 

Among all the decay modes, B c (ns) — > B c (ms) + it + tt in > m) is a favorable one for 
investigating the B c family because it is a strong-interaction process, moreover, it is also 
an ideal place to study the heavy hybrids \bcg > (\bc q >). Of course, the radiative decay 
B* — y B c + 7 is also a place to study the family 18(, but as our experimental colleagues 
suggest, at LHCb, the 7 detection would be difficult. Instead, the two charged pions are 
easy to be identified at LHCb detector. A rough estimate of the mass of B* in terms of the 
potential model indicates that the mass of B* is 6.36 GeV which is close to that given in 



Ref. (l9j. It only 80 MeV heavier than the ground state B c , thus the mode B* — > B c + it is 
forbidden by the final phase space. 

For the decay B c {ns) —> B c (ms)+7T+7T (n > m), the dominant mechanism is the t wo-g luon 



20|, the 



emission which eventually hadronize into two pions. As is indicated in the literature 
other minor mechanisms such as the subsequential pion emissions, may interfere with the 
amplitude induced by the two-gluon emission mechanism to change the lineshape of the 
differential width. But the total width is definitely determined by the two-gluon emission, so 
that our numerical results would give the decay width which can be measured in the future 
experiments at LHCb and Z-factory. The effects induced by the secondary mechanisms, may 
be measured at the Z-factory where a cleaner environment can provide an opportunity to 
conduct accurate measurements including the geometrical distribution of produced pions. 

In our calculations, we need to input several potential parameters to calculate the masses 
of the excited states of B c (ns) (n > 1). Unlike for the charmonia ip and T families, lack of 
data on their masses causes errors in our numerical results. 

We employ the QCD-multipole expansion (QCDME) method to calculate the correspond- 
ing widths. The specific potential forms for B c meson and hybrid beg are selected based on 
the lattice results. With the potential and concerned parameters we find that the mass of the 
ground hybrid beg is 7.63 ~ 7.64 GeV as the quark content be is in a spin-singlet. The widths 
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of B c (3s) — > B c (2s) + 7nr and B c (2s) — > B c (ls) + iT7r are not sensitive to the change of param- 
eters within the concerned range and it is similar to the cases for T(T(ns) — > T(ms) + tttt). 
While calculating the function we sum over the intermediate states of appropriate quan- 
tum numbers and truncate the expansion at K = 7. Our calculation suggests that the decay 
widths of B c (ns) — > B c (ms) + tt + tt can reach a few of tens of KeV which is of the same 
order as ip(ns) — > ip(ms) + tttt and T(ns) — > T(ms) + tttt, therefore with the luminosity of 
LHCb and the proposed Z-factory, there would be no problem to make relatively accurate 
measurements on such pion radiative decays. 

B c mesons were marginally produced at the LEP-I, as the luminosity of the proposed Z° 
can be at least three orders higher than that of LEP-I, there should be sufficient data on B c 
available. Moreover, there is a large phase space for the excited states B c (ns) (n > 1), so their 
production rates are not suppressed by the phase space and are similar to that for the ground 
B c , so that there should be a good chance to observe B c (ns) — > B c [ms) + tt + tt (n > m) at 
the Z° factory. As aforementioned, the background at the Z° factory is relatively small and 
a clean environment is expected. 

It is believed that LHCb and even TEVATRON possess sufficient database for observing 
such decays. In fact, B c production was first observed at TEVATRON, and with the energy 
and luminosity of LHCb, observation of B c (ns) — > B c (ms) + 7r + 7r (n > m) is definitely 
feasible. However, on other aspect, both LHC and TEVATRON are hadron colliders, so 
the background is much messier. With the efforts of the experts including theorists and 
experimentalists, it is already possible to clearly distinguish the signal from the background. 



Definite 
ground 



high quality generators are necessary for analyzing all possible sources of back- 



211 ] . By contrary, the background at the Z° factory is not so serious, i.e. the QCD 
contamination is relatively alleviated, even though it still exists. Detailed analysis on the 
possible background is a rather difficult job and usually is done by experts. When preparing 
the manuscript, we have consulted with our experimental colleagues about the possibility of 
observing such decays and analysis on the background, and here we can only make a very 
rough discussion. 

Therefore, for getting a better understanding of B c meson and its excited states, the Z° 
factory is definitely superior to the hadron colliders. 
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No doubt, the present work is still a primary effort to find the excited states of B c and 
study their structures, as well as that for the hybrid states \bcg >. One can be convinced 
that the order of magnitude of the numerical results is trustworthy, so that it is optimistic 
that measurements on such processes at LHCb and even the proposed Z-factory can be 
conducted. When the data are available, we will be able to further investigate the structure 
of the B c family and identify the mechanism(s) which governs the transitions. Then more 
precise theoretical works will be needed. 
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